Resilient Quantum Computation in Correlated Environments: 
A Quantum Phase Transition Perspective 
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We analyze the problem of a quantum computer in a correlated environment protected from 
decoherence by QEC using a perturbative renormalization group approach. The scaling equation 
obtained reflects the competition between the dimension of the computer and the scaling dimension 
of the correlations. For an irrelevant flow, the error probability is reduced to a stochastic form 
for long time and/or large number of qubits; thus, the traditional derivation of the threshold the- 
orem holds for these error models. In this way, the "threshold theorem" of quantum computing is 
rephrased as a dimensional criterion. 



A central result in the theory of quantum error correc- 
tion (QEC) is the "threshold theorem" Q. Even though 
QEC is a perturbative method [J, 0| , the "threshold the- 
orem" states that: Provided the noise strength is below a 
critical value, quantum information can be protected for 
arbitrarily long times. This remarkable theorem was first 
derived for stochastic error models Q or, more generally, 
for environments where the spatial and temporal corre- 
lations have an exponential decay Although there is 
still some controversy [5] , the theorem is widely accepted 
by a community where the foremost effort has been to 
extend it to correlated error models d, 0, B • 

The key difference between a correlated error model 
and a stochastic one is the absence of a local and time 
independent error probability. The main concept in ex- 
tending the theorem to a correlated environment has 
been the norm of an operator 0, 0, Q • A convenient norm 
of the interaction Hamiltonian can be used to bound the 
error probability and, eventually, to prove the thresh- 
old theorem. However, for some models the norm of 
the interacting Hamiltonian is too large to provide use- 
ful bounds 0. In this paper, we hence choose to follow 
different reasoning. 

We demonstrate that a large class of correlated error 
models is reduced to a simple stochastic error model in 
the asymptotic limit of large number of qubits or long 
time. Thus, in order to prove the resilience of the quan- 
tum information, we can fall back on the traditional 
derivation of the threshold theorem. Because the condi- 
tions for this fall back have clear parallels with the the- 
ory of quantum phase transitions [lol . 11|. we rephrase 
the threshold theorem as a dimensional criterion: (i) For 
systems above their "critical dimension" , the traditional 
proof of resilience is valid, and there are two regimes, or 
phases, as a function of the coupling with the environ- 
ment, (ii) However, when the system is below its "critical 
dimension" , the effects of correlations produce large cor- 
rections, and it is not possible to prove resilience by our 
arguments. 

In an elegant paper, Aharonov, Kitaev, and Preskill [8] 



(AKP) derived a new set of "threshold conditions" for a 
correlated noise model with instantaneous power law in- 
teractions between any two qubits of the computer. AKP 
proved that a version of the theorem holds if: (i) the cou- 
pling is below a critical value, and (ii) the interaction be- 
tween the qubits decays sufficiently fast when compared 
with the number of spatial dimensions of the computer. 
The first condition is very much what one expects for 
the "threshold theorem" . However, the second is very 
suggestive of our interpretation of the threshold theorem 
in terms of a quantum phase transition. Even though 
AKP considered an error model substantially different 
from ours (see below) , a particular case of our discussion 
reproduces AKP's conditions for resilience. 

The error model — We consider an environment de- 
scribed by a "non- interacting" field theory, i/o, allowing 
application of Wick's theorem. Furthermore, since an 
environment with a spectral gap would imply exponen- 
tially decaying correlations 0], we assume that the envi- 
ronment is gapless, has wave velocity v, and short time 
cutoff A"-'^. We focus on an interaction Hamiltonian that 
is local in the qubits. 



y = ^ Aq/q (x) CTq (x) , 



(1) 



where / is a function of the environment degrees of free- 
dom and fj are the Pauli matrices that parametrize the 
qubits. Finally, we consider (without loss of generality) 
the qubits to be arranged in a _D-dimensional lattice. 

The evolution of the system in the interaction picture 
during a QEC cycle is given by 



J7(A,A„) = Tte 



^'7o'^'itE.,„A»/»(x,t)<T„(x) 



(2) 



with A corresponding to the time of the syndrome ex- 
traction, Tt the time ordering operator, and /^(x, t) — 
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/.(x) 



At time A, the extraction of syn- 



dromes selects some terms of Eq. ([2|) as the evolution op- 
erator for that particular QEC cycle Hence, we 
naturally define a coarse-grained space-time grid of hy- 
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percubes, A x (vA) ' ^ (with z being the dynamical ex- 
ponent of the environment). 

The coarse grain grid is the frame upon which we de- 
velop our analysis. It provides the scale that separates 
two distinct noise regimes: intra- and inter-hypercube 
components. The key simplifying assumption is that in 
each volume (vA)^^^ there is only one qubit. In this 
case, the intra-hypercube part is simply the probabil- 
ity Ca of an error of type (Tq on that particular qubit. 
Two physical systems where this hypothesis should be 
immediately valid are acoustic phonons interacting with 
solid-state qubits and ohmic noise due to voltage/current 
fluctuations on qubits based on quantum dots or super- 
conducting devices. 

Our discussion is therefore divided into two parts. 
First, we demonstrate how to calculate this "stochastic" 
error probability, ea- Second, we evaluate how this error 
probability is changed by the inter-hypercube component 
of the noise. 



Defining Eq. — Consider, for instance, that we know 
from the syndrome that a particular qubit, say Xi, suf- 
fered an error (Tq, in the time interval labeled 0. To lowest 
order in Xa, the QEC code disentangles the qubit from 
the environment. Thus, 

Ua (Xi, Aq) « -iAa / dt/a(Xi,t), (3) 

Jo 

is the evolution operator associated with that particular 
qubit, where we used that ctq, (A) (t) = 1 [13|. 

Of course, lowest order perturbation theory is not al- 
ways justified. Hence, before QEC can be argued to be 
effective, it is important to estimate how higher order 
terms change Eq. ([3]) . A very direct approach is to use 
the perturbative renormalization group (RG). 

To derive the RG equations, we supplement ^ with 
the next higher-order terms allowed by the QEC code 



{)a(xi,Aa) « -iAq [ /a (xi,i) - i |eQ/37 1 A/3 A^ o-Q (A) Tt f d^i (it2 //3 (xi , ii) (xi, ta) cr^j (ti) ct^ (t2) 
JO ^ Jo 

+ - AaA^CTa (A) Tt / dilrfi2rfi3/a(xi,tl)//3(xi,<2)//3(xi,i3)o-Q(^l)cr/3(i2)cr/3(t3), (4) 
Do Jo 
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where Eq/j^ is the antisymmetric tensor. Since QEC com- 
pletely removes the qubit variable from the problem, we 
can consider Eq. ^ as a field theory problem in itself. 
It is therefore straightforward to derive the lowest order 
terms to the /3- function of ■Ua(xi), 



^ = 5/37 (^) ^/3-^7 + W AaA^, 



(5) 



where g and h are functions specific to a particular envi- 
ronment and d£ = dA/A JJJ. Note that the integration 
of Eq. ([5]) is equivalent to summing an infinite series of 
terms in the perturbative expansion for Va (xi ) . Because 
there is only one qubit in the hypercube, the problem of 
calculating Eq, involves only /(xi) and so was reduced to 
an impurity problem. 

To further proceed with the argument, we must ask 
that the renormalized value of Aq, at the (wA)^^scale, 
A* , be a small parameter. In that case, it is appropriate 
to use the evolution Eq. ([3]) with Aq replaced by A* . 

The conditional probability of having an error of type 
a in a particular hypercube labeled by xi and t — has 
the general form [12] 

P(...;a,xi; ...) w (...-u^ (xi. A* ) ...£>« (xi. A* ) ...) . 

Hence, following the discussion in Ref. we define the 



operator that gives the probability of an error as 

(xi, a;) « vi (xi, a;) Va (xi, a;) . (6) 

We can now readily separate the effects of correlations 
into their intra- and inter-hypercube parts. Using Wick's 
theorem, we obtain 

i;2(xi,A;)«eQ + (A:Af : I/„ (xi, 0)|' : , (7) 

where the intra-hypercube part is 

/•A i-A 

e^^iXlf dt, dt2 (/I(xi,ii)^(xi,i2)) (8) 
Jo Jo 

and :: stands for normal ordering. An important remark 
is that Ca is from this point on a numerical factor set by 
A* at the scale {vAy^. 

For later convenience, we re-write Eq. ([7]) in a slightly 
different form. Defining the operators 



^^a(xi,0) 



(a;a)= 



l/a(xi,0)r 



(9) 



we re- write the operator for the probability of an error as 
the product v^. (xi, A* ) = [1 + Fa (xi, 0)]. By direct 
calculation, or simply by the unitarity of the probabil- 
ity, it is also straightforward to write the operator for 
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the probability of not having an error as Vq (xi. A*) = 
[1 - J2a ^a] -^0 (xi, 0), where 

^^(a^a)' :|/^(xi,0)|': 
Fo(xi,0) = 1- ■ (10) 

^ L^I3=x,y,z^l3 

Probability function and scaling — Now that we sepa- 



rated the probabihty into an intra- and inter-hypercube 
component, we can seek to write the probabihty of a 
computer evolution with m errors after N QEC cycles. 
A particular case is helpful in understanding how to pro- 
ceed. Suppose we want to calculate the probability that 
R qubits suffer m <C NR errors of type a. Using Eqs. 
(l7])-(l9l), it is straightforward to write 



P" — 7) 
^ m — Prn 



{vAf/' (vA)^^" Jo 



ir ^ ( [n ^o(xc, tc)] [1 + F„(xi, ti)] ... [1 + ^^.(xn., t™)] 



where we integrated over all possible grid positions, 
{'Xj,tj), C denotes the set of remaining hypercubes, and 

Pm = Ea)^^"™ (Ca)'" [lB| . 

We now organize the expectation value of Eq. (fTTjl in 
powers of (A*, ea) and invoke Wick's theorem again. The 
first term is just the stochastic contribution to the prob- 
ability, 

(12) 

The next term is typically of the form 

J k=l (^^) 

Thus, the fundamental role of the scaling dimension of 
Fa now becomes clear. If dim [/„] = 5a, then Eq. ^ 
implies that dim \_Fa\ = 25a- Since Hq is assumed to be 
non-interacting, the two point correlation function has 
the general form 

{Fa (Xj, t^) Fa (Xj, tj^) ^ J~ ^li-i Xj I , tj\ ^ ^ 

For large N or i?, we can study the stability of the ex- 
pansion of Eq. ([TTl) in powers of A* using the traditional 
scaling theory [lO|. The simplest way to proceed is to 
apply the transformation Xjj &X;j and ti^j b^t'i.j^ 
with b = e''^, to Eq. (fT3|) . It is then straightforward to 
obtain the scaling equation for A* , 

d\* 

_^ ^ iD + z-dim[Fa])K. (14) 

A similar argument can be applied to all the remain- 
ing terms of Eq. (fTTjl . The critical condition is when 
D + z = 25. In this particular case, correlations be- 
tween hypercubes introduce logarithmic corrections to 
the stochastic part of Eq. pT|) . 



An irrelevant flow for A* indicates that the system 
is above its "critical dimension". Since correlations be- 
tween hypercubes produce small corrections in compari- 
son to Eq. (jl2p , the probability distribution at long time 
and/or large number of qubits has essentially a stochas- 
tic form. Corrections to this form can be systematically 
calculated by perturbative expansion in A*. Thus, we 
can fall back on the usual proof of resilience for quantum 
information 01 : There will be two "phases" separated 
by a A"'*'''''': (i) For Aq < A"'"''''', the information can 
be protected for arbitrarily long times. Hence, the com- 
puter and the environment are disentangled by the QEC 
code, (ii) In contrast, for A^ > A"''*"^''\ the computer 
and the environment can not be disentangled and deco- 
herence will take place. 

This scenario strongly resembles the theory of quantum 
phase transitions. There are several ways to pursue this 
analogy; here we present two. First, we just showed that 
for systems above their critical dimension, the second 
term of the r.h.s. of Eq. ([7]) is unimportant at asymptot- 
ically large scales. In other words, when we calculate the 
probability P (....; a , Xi ;...), it is a good approximation to 
replace the operator Va (xi. A* ) by its "mean field" (per- 
turbative) value Cq. The bare value of Aq is the only pa- 
rameter that determines whether QEC can succeed and 
thus plays a role analogous to temperature. This is pre- 
cisely what one expects for systems above their upper 
critical dimension, where the transition has a mean field 
character [lo| . 

Second, an explicit analogy with a statistical mechan- 
ical problem is also possible. Eq. (fTT|) can be thought 
of as the partition function of a gas of fictitious parti- 
cles, where A* is the particle fugacity. Fa creates a par- 
ticle, and Fq introduces the vacuum fluctuations. An 
example is quite illuminating: Consider the simple case 
of a bosonic bath with (0, tt) representing the conju- 
gate fields. For a Z) = 1 computer and a noise model 
V = X J2j cos [5 ■ (j) (j)] (t|, Eq. (fTTj) is precisely the parti- 
tion function of a 2-dimensional Coulomb gas. This case. 
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then, has a quantum phase transition in the KosterHtz- 
Thouless universahty class as a function of S. 

On the other hand, going back to Eq. (|14p . a relevant 
flow for A* indicates that inter-hypercube correlations 
produce contributions to the probability function that 
scale in the same fashion as the intra- hypercube part. In 
this sense, there is no "de facto" separation of scales or, in 
other words, the problem is inherently non-perturbative 
in the coupling with the environment. It is therefore not 
possible to prove resilience by our arguments. 

Changing dim [Fa] — We showed previously that 
small changes in the QEC code can dramatically reduce 
the effects of correlations between the hypercubes. By 
simply applying logical NOTs and phase-NOTs in each 
QEC cycle, the dimension of Fa becomes larger at the 
cost of increasing Eq. To understand this fact, revisit 
Eq. ([3]); for illustration, consider a = z and apply a logi- 
cal NOT half way through the cycle: 



(xi. A) w zAz 



dt 



A/2 



A/2 



dt 



/.(xi,i). 



Following the same steps as before, we obtain that 
F^ (xi,0) oc : \dtfz (xi,0) \^ :. If n logical NOTs and 
phase-NOTs are performed, we find that dim [Fa] — 
2 {da + nz). Therefore, for a given noise model, one can 
always engineer an irrelevant flow. 

Connection to the AKP results — In Ref. Q a family 
of long-ranged noise models with interactions between 
qubits was studied. AKP considered a power law inter- 
action between any two qubits at positions xi and X2 of 
the computer with strength Fjxi — X2|^^'^. Clearly, one 
could start from their noise model and use a Hubbard- 
Stratonovich transformation to arrive at ours. The re- 
verse is also true: starting from our model, one could 
integrate out the environment field and arrive at the ef- 
fective interaction between qubits that AKP considered. 

Nevertheless, the two papers deal with opposite limits 
of this same model. The crucial difference is the wave ve- 
locity of the environment, v. AKP considered the limit of 
instantaneous interactions, — > oo. This implies that all 
the qubits of the computer are contained in the volume 
(vA)^^^, in contrast to our assumption that one qubit 
is in that volume. Interestingly, the AKP limit leads to 
an important simplification when combined with QEC: 
Since errors are detected by QEC, there are no mem- 
ory effects due to correlations between qubits at different 
QEC steps. Hence, their error model corresponds in our 
analysis to a dynamical exponent z = 0. With this, our 
criterion for the possibility of resiliency is exactly the 
same as AKP's, even though the problem and methods 
used are inherently different. We speculate, therefore, 
that there is a more general scheme that encompasses 
both papers. 

Remarks and conclusions — In hindsight, the results 
that we obtain are dramatically clear and simple. QEC 



is a perturbative method. Therefore, the "threshold con- 
dition" should be a statement about when a perturba- 
tive analysis is valid. In field theories this is a relatively 
straightforward question. However, in a quantum com- 
puter this is not so obvious due to the presence of the 
qubits. The remarkable result of QEC is that some of 
the quantities relevant for the threshold condition depend 
exclusively on the environment variables. It is, then, pos- 
sible to derive criteria based only on the field theory that 
describes the environment. 

Our discussion brings to light an interesting parallel 
with the theory of quantum phase transitions. In fact, 
we derived a condition that strongly resembles the def- 
inition of the "upper critical dimension" of a quantum 
problem. When the system is above this dimension, the 
usual "mean field" derivation of the threshold theorem is 
applicable. Hence, as a function of the bare coupling with 
the environment, there is a transition from the states of 
the qubits and the environment being disentangled to 
their being entangled. 

An important question that remains open is if there 
is also a "lower critical dimension" , namely a criterion 
for the impossibility of proving resilience. If the lower 
and upper critical dimensions are not the same, the in- 
termediate cases would require a proof of resilience sub- 
stantially different from the "mean field" approach used 
here. An example of such an approach is the use of the 
norm of the interaction Hamiltonian by AKP. However, 
the fact that they found the same "critical dimension" 
as we have found suggests that the "upper" and "lower" 
critical dimensions may coincide. 
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